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ABSTRACT 

We study the radio/X-ray correlation in Cyg X-3. It has been known that the soft and hard 
X-ray fluxes in the hard spectral state are correlated positively and negatively, respectively, 
with the radio flux. We show that this implies that the observed ~1—100 keV flux (which is a 
fair approximation to the bolometric flux) is completely uncorrelated with the radio flux. We 
can recover a positive correlation (seen in other sources and expected theoretically) if the soft 
X-rays are strongly absorbed by a local medium. Then, however, the intrinsic X-ray spectrum 
of Cyg X-3 in its hard state becomes relatively soft, similar to that of an intermediate spectral 
state of black-hole binaries, but not to their true hard state. We also find the radio spectra in 
the hard state of Cyg X-3 are hard on average, and the flux distributions of the radio emission 
and soft X-rays can be described by sums of two log-normal functions. We compare Cyg X-3 
with other X-ray binaries using colour-colour, colour-Eddington ratio and Eddington ratio- 
radio flux diagrams. We find Cyg X-3 to be spectrally most similar to GRS 1915+105, except 
that Cyg X-3 is substantially more radio loud, which appears to be due to its jet emission 
enhanced by interaction with the powerful stellar wind from the Wolf-Rayet donor. 

Key words: radiation mechanisms: non-thermal - radio continuum: stars - stars: individual: 
Cyg X-3 - stars: winds, outflows - X-rays: binaries. 


1 INTRODUCTION 


Cyg X-3 is an binary system possessing a number of unique and 
puzzling cha racteristics. It is the brigh test radio source among X- 
ray binaries dMcCollough et alJll999l), showing extremely strong 
radio outbursts and resolved jets fe.g.. lMartt. Paredes & Peracaulal 
l200ll : Mk>duszewski et al J|200 11) . It is the only X-ray binary that is 
certainly powered by accretion for which high-energy y-ray emis¬ 
sion has been detected at a high statistical significance (Tavani et al. 
120091 : iFermi LAT Collaboration! [20091) . It is the only known bi¬ 
nary in the Galaxy containing both a compact object and a Wolf- 
Rayet star 1 vanKeAwijketal. 1 19921 . 1 19961 : Ivan Kerkwiikl 1 19931 : 
iFender. Hanson & Poo lev! 1 1 ~ ). Its orbital period of P ^ 0.2 d 
is unusually short for a high-mass binary, which indicates a past 
spiral-in episode during a common-envelope evolutionary stage. 
Given the value of P and estimates of the masses and mass-loss 
rate, the compact object orbits within the donor photosphere. 

In spite of its discovery as early as in 1966 (iGiacconi et alj 
119671) . the nature of its compact object remains uncertain, due 
to the lack of fully reliable determination of the mass functions 
and inclination. Based on the currently available information, 
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IZdziarski. Mikolaiewska & Belczvnskl d2013l) found the range of 
the compact-object mass of 2.4 /jJMq, allowing either a neutron 
star (NS) or a low-mass black hole (BH). The nature of the compact 
object is of significant interest, in particular because Cyg X-3 is a 
certain progenitor of a doubl e compact system, after its Wolf-Rayet 
star explodes as a supernova l lBelczvhski et al.l2013l) . The presence 
of a BH o ver an NS is favoured by considering the X-ray and radio 
emission iHjalmarsdotter et al. ) l2008j. I2OO9I: ISzostek & ZdziarskH 


20081 ISzostek. Zdziarski & McColloughl 200a. herea fter SZM08: 
Kolio nen et al. 2010h . Zdziarski. Misra & Gierlinskil j 20 1 ()t) have 


shown that the differences between the shape of the X-ray spectra 
of Cyg X-3 in its hard spectral state from those of confirmed BH bi¬ 
naries could be accounted for by Compton scattering in the strong 
stellar wind from the donor. That model als o co uld acco unt for 
the lack of high frequencies dAxelsson. Larsson & HialmarsdotteiH 
120091) in the power spectra of Cyg X-3. 


The distance to Cyg X-3 has been estimated as D ^ 7-9 kpc 
dLing, Zhang & Tang] 1 20091 : iDickevI 1 19831 : IPredehl et alj|2000h : at 
this D, its observed (i.e., not corrected for intrinsic absorption) 
bolometric X-ray luminosity is ~ 10 38 e rg s -1 . Its X-ray spec - 
tra have been classified into six states by iKolionen et alj J2010h . 
based on Rossi X-ray Timing Explorer ( RXTE ) pointed observa¬ 
tions. Here, we adopt a simpler classification into three states, hard, 
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Figure 1 . (a) The correlation between the daily-averaged energy fluxes in the 3-5 keV (ASM) and 14-50 keV (BAT) ranges. The hard and soft states are 
defined here by F( 3-5 keV) < 0.5 keV cm' 2 s 1 and f { 14 50 keV) < 1.2 keV cm -2 s 1 . respectively. These boundaries are marked by the dotted and dashed 
lines, respectively. An intermediate state corresponds to both the 3-5 keV and 14-50 keV fluxes above the respective boundaries. Only points with statistical 
significance > 3cr (for each of the fluxes) are shown, (b) The same except that the 2—4 keV flux from MAXI is used. The adopted hard-state condition (the 
dotted line) is here F( 2-4 keV) < 0.4 keV cm' 2 s' 1 . 


intermediate and soft, based on the relationship between the 3-5 
keV and 14-50 keV fluxes. An unresolved issue is the value of the 
absorbing column wi t hin th e stellar wind of the donor. For example, 
iHialmarsdotter et aD d2008i ) found that the hard-state spectral data 
from INTEGRAL can be fitted with either a high or low absorbing 
column. The two resulting models then significantly differ in both 
the spectral slope of the soft X-rays and in the bolometric flux. In 
the former and latter models, the intrinsic spectrum is soft, with the 
photon index of F >2, and hard, with T < 2 (defined by the pho¬ 
ton spectrum, cc £' r ). Only the latter behaviour corresponds to_the 
standard definition of the hard state, e.g. jRemillard & McClintockl 

d2006h . 


Correlated radio X-ray behaviour of Cyg X-3 was discov¬ 
ered in the soft (1-6 keV) and hard (20-100 keV) X-rays by 
IWatanabe et al.l d 19941) and lMcCollough et alj j 1999^ respectively. 
The c o rrelations were then studie d by iGallo. Fender & P oolev 
d2003l) . iHialmarsdotter et akl ( 2008 > and SZM08. Curiously, the 
soft X-rays are positively correlated with the radio flux in the 
hard state whereas the hard X-rays are anticorrelated. This dif¬ 
fers, e.g., from another accreting high-mass binary, Cyg X-l, in 
which both the soft and hard X-rays are positively_correlated with 
the radio emission in the hard state dZdziarski et al.ll201 ll : here¬ 
after ZSPL11). This raises the issue of how the bolometric flux 
(dominated by the X-rays, emitted most likely by the accretion 
flow) is correlated with the radio emission, originating in the jets 
of this system. Theoretically, a positive co rrelation is expected 
jHeinz & Sun vaevf 20031 : Merloni et alj2003l ) in a steady accretion- 
flow-jet system. Generally, this is because both the accretion lumi¬ 
nosity and the jet power should be positively correlated with the 
accretion rate, and the radio flux is positively correlated with the 
jet power faeinz & Sunvaevi2003h . 


In this work, we study in detail the relationship between the 
radio emission and X-rays. In particular, we estimate the bolomet¬ 
ric flux of the system, its relation to the radio flux, and the radio 
spectra. We then consider implications of our findings for the na¬ 
ture of this accreting source, and compare Cyg X-3 to other X-ray 
binaries. 


2 THE RADIO AND X-RAY MONITORING DATA 

We study here radio monitoring data at 15 GHz from the Ryle Tele¬ 
scope, which cover MJD 49231-53905, and the Arcminute Mi- 
crokelvin Imager (AMI), MID 54612-57211. The combined data 
set contains 76776 measurements. The AMI Large Array is the 
re-built and reconfigured Ryle Telescope. IPoolev & Fended <1 1 9971) 
describe the normal operating mode for the Ryle telescope in the 
monitoring observations; the observing scheme for the AMI Large 
Array is very similar. The new correlator has a useful bandwidth of 
about 4 GHz (compared with 0.35 GHz for the Ryle), but the effec¬ 
tive centre frequency is similar. The data are subject to variations 
in the flux calibration of $ 10 per cent from one day to another; 
we hereafter assume a 10 per cent fractional error. Also^we use 
the 2.25 GHz and 8.3 GHz monitoring data dWaltman et al .11 1 99fil 
Il995l[l99^) from the Green Bank Interferometer (GBlfl The data 
are for MID 50409-51823 and contain 11022 measurements. 

We average all of the data used here over 1 day (except for 
the BATSE, see below). We have found that 1 day is a time scale 
which both strongly reduces the measurement errors and averages 
over the orbital modulation (with the orbital period of 1/5 d), see 
also IZdziarski et akl d2012l) . At the same time, it retains most of the 
correlated radio and X-ray variability. 

For soft and medium X-rays, we use daily-averaged 
monitoring data from the All-Sky Monito r (ASM; 
iBradt. Rothschild & Swanki i 19931 : L evine et al. 19%!) on board 
RXTE, giving count rates in three channels at energies of 1.3-3 
keV, 3-5 keV and 5-12.2 keV, spanning MID 50087-55915. We 
also use the data from the Monitor of All-sky X-ray Image (MAXI; 
Matsuoka et al. 120091) on board of International Space Station 
(ISS). We use daily averaged data during MID 55058-57401. The 
datc0 are given as photon fluxes in three channels, 2-4, 4-10 and 
10-20 keV. These are converted from the count rates by scaling to 
Crab assuming its photon spectrum of the form of 


1 ftp://ftp.gb.nrao.edu/pub/fghigo/gbidata/gdata/ 

2 http://maxi.riken.jp/top/ 
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F(E)/E = A exp(-N H cr H )(E /1 keV)' 


(1) 


where the H column density is Ah = 3.5 x 10 21 cm' 2 , <x H is the 
bound-free cross section per H atom averaged over the cosmic com¬ 
position (calculated using wabs in xspec. lAmau 3 l 19961) . A = 10 
cm' 2 s' 1 , and T = 2.10. We convert the MAXI data to the energy 
fluxes using the above. For other data sets used here, we use the 
above form of the Crab spectrum except for taking into account a 
break to F = 2.24 at > 100 keV. 

For hard X-rays, we use daily-averaged data from the 
Swift Burst Al ert Telescope (B AT; iBarthelmv et alJ 120051 : 
iMarkwardt et al.| [2005: Kri mm et~aill20 1 3l) survey data. They were 
retrieved from the HEASARC public archive and processed using 
the BAT_IMAGER software (Segreto et alj]20ld) . They form a 14— 
195 keV 8-channel count-rate light curve, similar to that used in 
Izdziarski et alJ I 2 OI 2 I) . but covering now a substantially longer in¬ 
terval of MJD 53348-57279. The energy fluxes for the ASM and 
BAT data have been obtained by using scaling to the Crab spec- 
trum as given above. We note that lin’t Zand. Jonker & Markwardti 
d2007h have estimated that, if the photon index of the actual source 
differs from the Crab index by 1 or the column density is up to 10 
times larger, the ASM energy-range flux differs by at most 30 per 
cent. We thus expect similarly minor effects on the scaling of the 
MAXI and BAT data. 

We also use data from the Burst and Transient Source 
Experiment (BATSE; llTarmon et al.l 120021) on board of Comp¬ 
ton Gamma Ray Observatory ( CGRO ). These monitoring data 
have been _obtained using the Earth occultation analysis tech- 
nique (Harmon et al. 2002). We use the same data as used by 
iFIialmarsdotter et alJ 1 2008 ) and SZM08, which are given as 3- 
day average energy fluxes in the 20-50, 50-100 and 100-230 keV 
bands for the interval of MJD 48363-51686. 

In Appendix [A] we present and discuss the histograms of the 
flux distributions and the light curves for the X-ray monitors and 
for the Ryle/AMI. 


3 THE RADIO/X-RAY CORRELATION 

Fig- H defines the spectral states of Cyg X-3 used in this paper 
via the relationship between the soft and hard X-rays, similarly 
to Izdziarski et alJ 112012l) . The horizontal dashed line in both pan¬ 
els delineates the upper boundary of the soft state, with F(14— 
50)keV) <1.2 keV cm' 2 s' 1 . The vertical dotted lines delineate 
the upper boundary of the hard state, with F(3-5) keV) < 0.5 keV 
cm' 2 s' 1 for the ASM data, and F(2-4) keV) < 0.4 keV cm' 2 s' 1 
for the MAXI data. The intermediate state is defined for the fluxes 
above both boundaries. 

Fig.[2]shows the energy dependence of the 15-GHz/X-ray re¬ 
lationship based on the ASM/BAT (top) and MAXI/BAT (bottom) 
data. Each data point corresponds to a daily average of the X-ray 
flux and the corresponding average radio flux. In Figs. 0a-c) and 
(e-g), the data points showing a strong positive correlation seen on 
the left correspond to the hard state. The dashed lines in Figs.[2jb,e) 
show the upper flux limits of the hard state in the 3-5 and 2^1 keV 
bands, respectively, defined above. 

As found by SZM08, as the soft X-ray flux increases above 
the upper hard-state limit, the radio flux becomes highly variable, 
dropping to very low values and increasing to very high values, cor¬ 
responding to major radio flares. Fig.[2]does not show the temporal 
behaviour (shown in SZM08; see also Appendix[A]l, but we do see 
the radio flux varying by 4 orders of magnitude while the soft X-ray 
flux changes only within a factor of ~2. 


As we also see in Fig. [2] the slope of the radio/X-ray corre¬ 
lation changes with the energy, becoming increasingly steep with 
increasing energy up to ~10 keV, and then changing its sign to 
negative in the 14-50 and 20-75 keV bands. The vertical dashed 
lines in Figs.[2jd,h) show the lowest flux in those bands correspond¬ 
ing to the hard state. The change of the sign does happen between 
the ASM and BAT energy ranges, as we have checked by looking 
at the 14-20 keV band, which already shows a negative correla¬ 
tion (though with more scatter). Since that dependence is similar to 
those in the 14-50 and 20-75 keV bands, we do not show it here. 

We have found the steepening of the correlation with the en¬ 
ergy in the hard state in the simultaneous ASM/GBI data, at radio 
frequencies of 2.25 and 8.3 GHz, to be very similar to those shown 
in Figs. [2]A-c,d-g). Also, we have correlated the BATSE and 15 
GHz data. We have found an anticorrelation at 20-100 keV range 
in the hard state similar to those shown in Figs. [2jd,h), but with 
much higher statistical scatter. The 100-230 keV BATSE data have 
then very limited statistics, given the spectrum of Cyg X-3 being 
strongly cut off at those energies. Given the similarities of those 
results to those using the ASM, BAT and MAXI and the 15-GHz 
data, we do not show them here. 

The existence of the anticorrelation between radio and hard 
X-rays in the hard state raises the issue of how the bolometric lumi¬ 
nosity is correlated with the radio flux. The observed EF(E ) spectra 
in the hard state peak at ~20 keV (e.g., SZM08), and the spectral 
slope in soft X-rays is F ~ 1, which gives rise to the classification 
of this state as ‘hard’. Thus, the anticorrelation takes place in the 
spectral region in which at least ~ 1/2 of the observed bolometric 
flux is emitted. 

Figs.[3ja ,c) show the relationship between the 1.3-100 and 2- 
100 keV flux, respectively, and the radio flux. As it can be seen 
in the spectra of, e.g., SZM08, these X-ray energy ranges contain 
most (though not all) of the observed X-ray energy and thus they 
give fair approximations to the bolometric flux. (For the ASM/BAT 
data, we have neglected here a minor contribution from the 12.2-14 
keV band, not covered by either ASM or BAT.) We see that there 
is hardly any correlation, with the data just forming a long vertical 
strip. 

Obviously, there is substantial absorption in Cyg X-3, which 
need s to be taken into account. Its amount is, however, rather uncer¬ 
tain. iHialmarsdotter et all j200 8f) found that their hard-state data, 
from INTEGRAL , can be fitted by either weak absorption, with the 
resulting in trinsic slope below 10 keV of T < 2, or strong one, 
with T > 2. Szostek & Zdziarskil l2008i) found that their hard-state 
BeppoSAX data can be fitted with a weak absorption (T < 2). As we 
find here, only relatively strong absorption, corresponding to T ^ 2, 
gives reasonable indices of the hard-state correlation of the intrin¬ 
sic bolometric flux, Fb 0 i, with the radio flux, Fr, Fr oc FT. Thus, 
the hard state of Cyg X-3 would then be similar to an intermediate 
state of BH binaries. 

The correlation index for most BH low-mass X-ray binaries 
(LMXB s) in the hard state for the 3-9 keV X-ray range is p ^ 
0.6-0 .8 dCorbel et al.ll2003. 120031. Eooi I2008L 12013 d: iGallo et aD 
120031) . Then, some LMXBs have bee n found to have p - 1.4, in 
particular the BH LMXB H1743-322 JCoriat et alj|201 ll) . with the 
two branches intersectin g at L(l-lOkeV) ^ 10 38 erg s' 1 , s ee fig. 
9 in lCorbel et al.l (20 1 3al ). As discussed in lCoriat et al] (201 ll) , the 
two branches of the correlation correspond to the jet being powered 
by either radiatively inefficient {p — 0.7) or radiatively efficient 
(p - 1.4) accretion. 

Figs. [3jb,d) show the correlation between the radio flux and 
the intrinsic, unabsorbed, broad-band X-ray flux for a choice of the 


© 2016 RAS, MNRAS flOfl.lJUTTl 




























































4 A. A. Zdziarski et al. 


—3 

£ 


N 

X 

o 

10 


0.05 0.1 0.5 1 0.1 0.5 1 0.5 1 5 0.05 0.1 0.5 1 


F(1.3-3 keV) 


F(3-5 keV) F(5-12.2 keV) F( 14-50 keV) 



i i i |-1-1-1—i—i i i i |- p—i i i i |-1-1- —i—i—i i i i |-1-n -1-1—i—| i i i i |-1-1—i—| i i i i |-r 



0.1 1 1 1 0.05 0.1 0.5 1 


F(2-4 keV) F(4-10 keV) F( 10-20 keV) F(20-75 keV) 


Figure 2. The relationship between the 15 GHz flux and the X-ray flux in four energy bands based on (top panels) the ASM and BAT data and (bottom panels) 
the MAXI and BAT data, in units of keV cm -2 s _I . Only X-ray fluxes with significance > 3 cr are plotted. The dashed lines in (b) and (e) show the approximate 
upper boundary of the soft X-ray flux of the hard state, above which there are both the intermediate and soft states, see Fig. |T] The dashed lines in (d) and (h) 
show the approximate lower boundaries of the hard X-ray flux in the hard state, at 1.2 and 0.8 keV cm -2 s -1 , respectively. The vertical extent of these lines 
is the same as in (b) and (e), and it corresponds to the range of the radio fluxes at the boundary of the hard state. We see the sign of the hard-state correlation 
changing around ~ 10-20 keV. The strong radio flares occur only after the hard X-ray flux drops below its hard-state boundary (the dashed line), but then the 
radio flux in the soft state becomes (weakly) positively correlated with the hard X-rays. 
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Figure 3. Left panels (a-b): the relationship between the 15 GHz flux and the 1.3-100 keV flux, which is a good approximation to the bolometric flux, (a) The 
bolometric flux based on the ASM and BAT data as observed, (b) The flux corrected for absorption. Only X-ray fluxes with significance > 5<r are plotted.The 
dotted line shows the power law fit to the hard state. Right panels (c-d): the same as (a-b) but based on the MAXI and BAT data. The vertical dashed lines in 
(b) and (d) show the assumed upper boundary of the hard state at 15 and 12 keV cm -2 s -1 , respectively. See Section[3]for details. 







Figure 4. Left panels (a-b): the average spectra based on the ASM and BAT data in the hard (red triangles), intermediate (green squares), and soft (blue circles) 
states; (a) observed, and (b) corrected for absorption. Left panels (c-d): the same but based on the MAXI and BAT data. See Section[3]for details. 


intrinsic absorber yielding p ~ 1.5. We use the ionized absorber 
model zxipcf of lReeves et alj (120081) in xspe c, which is based on 
the tr eatment of ionized absorption by xstar (Bautista & Kallman 
1200 ll) . It is characterized by a column density, N, and an ionization 
parameter, £ = Lx/inr 2 ), where n is the medium density and r is 
the distance from the ionizing source to the medium. We assumed 
N = 10 24 cnT 2 and £ = 10 28 e rg cm s -1 , which values are similar 
to those fitted in tzdziarski et all i2012h to the average RXTE X-ray 
spectra for the orbital phases around the superior and inferior con¬ 
junctions. In addition, we took into account interstellar absorption 


with Ah = 1.5 x 10 22 cm 2 lLauaue et al.lll972l : IChu & Biegind 


119731: iDickev & Lockmanll 1990K 


We have then created average X-ray spectra in the three 
spectral states, defined as in Fig. |T] based on the simultaneous 
ASM/BAT and MAXI/BAT data. The spectra as observed and the 
intrinsic ones, i.e., corrected for absorption, are shown in Fig. [4] 
We see that the intrinsic spectra in the hard state are almost flat in 
EF(E) below ~10 keV, with T ~ 2. 


The slope of the observed radio/X-ray correlation correspond¬ 
ing to the hard-state ASM/BAT spectrum shown in Fig. |4jb), and 
fitted as described in ZSPL11, is p ^ 1.5 ± 0.1, which we show by 
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the dotted line in Fig. [31b). Given the overlap of the points in differ¬ 
ent states, the value of p also bears a substantial systematic error. 
This slope is very similar to that found in the lower-branch (radia- 
tively efficient) radio/X-ray correlation of LMXBs, see above. The 
corresponding slope for the MAXI/BAT data, shown in Fig. [3jd), 
is somewhat higher, which appears to be due to the lower boundary 
of the MAXI data being 2 keV compared to about 1.3 keV for the 
ASM data and the steepening of the correlation with the increasing 
energy up to ~10 keV, found here. 


4 RADIO SPECTRA 

In order to get more insight into the nature of the hard state, we 
consider the radio spectra based on three measurements, at 2.25, 
8.3 and 15 GHz. We calculate the spectral index, a, between two 
frequencies, defined by F v oc y“. We show the results in Fig. [5] We 
identify the hard state using the simultaneous ASM data, by our 
standard criterion of F(3-5 keV) < 0.5 keV cm -2 s _1 , which we 
show by the black symbols. Note that some points not identified in 
this way due to the lack of simultaneous ASM data also belong to 
that state. We see that the hard state corresponds to the range of the 
2.25-GHz fluxes of =35—150 mjy and a ~ -0.1 to +0.6, with the 
average of (a) — 0.2. Thes e radio spectra are harder than those of 
Cyg X-l, with a ^ 0 iFender et al.l200d) . However, they are similar 
to those of those exhibited by some BH binaries. For example, the 
LMXB GX 339^1 has showna_in the hard state during the 2010-11 
outburst saturated at = 0.4 ( ICorbel et alj2013b , see their fig. 2). 

Thus, the hard state in Cyg X-3 is characterized by hard ra¬ 
dio spectra, similar to the case of BH binaries. Such spectra are 
emitte d by jets with the partiall y synchrotron self-absorbed emis¬ 
sion jBlandford & Konlglll979l) . On the other hand, the soft state 
shows a wide range of a, with majority of them with a between 
-1 and 0. The emission with a < 0 is likely optically thin syn¬ 
chrotron. However, the presence of spectra with a > 0 during the 
intermediate and soft states, see Fig. [5] indicates that synchrotron 
self-absorption is also taking place in those states. 


5 DISCUSSION 

5.1 The intrinsic radio/X-ray correlation index 


We have found, based on the ASM, MAXI, BAT and Ryle/AMI 
data, that the Cyg X-3 radio flux in the hard state is highly vari¬ 
able and completely uncorrelated with the measured bolometric 
flux, as shown in Figs. [3ja,c). However, a positive correlation 
is recovered when taking into account absorption of soft X-rays 
provided it is strong enough, as shown in Figs. [3jb,d). We have 
found that an absorber similar to that fitted to the RXTE pointed 
data in lzdziarski et alj J2OI2I) gives a reasonable slope of the Fb 0 i- 
F(15 GHz) correlation, p ^ 1.5, similar to those found in the kiwer 
branch of the LXMB radio/X-ray correlation ( Corbel et al .11201 3 hi ) 
and in Cyg X-l in the hard state jZdziars ki il2012l) . 

However, this index can be reduced due to an effect_follow- 
ing from the physics of a cont inuous jet (e.g. jBlandford &~ Konigl 
ll979L [Heinz & Sunvaevll2003l) coupled with the free-free absorp¬ 
tion in the stellar wind. Namely, the height above the binary plane 
of a jet at which the bulk of partially synchrotron self-absorbed 
emission at a given v emerges increases with the increasing F v . 
This owes to the associated increase of the density of the jet rel¬ 
ativistic electrons, which in turn increases the synchrotron self¬ 
absorption optical depth within the jet. For the standard model of 



Figure 5. The spectral index (F v cc v a ) based on two radio frequencies 
vs. the flux of the lower frequency for 1-d average GBI and Ryle data, (a) 
2.25-8.3 GHz; (b) 2.25-15 GHz. The hard state points identified by the 
simultaneous ASM data is shown by black symbols. Only the 2.25. 8.3 GHz 
and 3-5 keV flux measurements with the significance of > 3 cr have been 
taken into account. 


iBlandford & Konigll Jl979t) and assuming that the radio emission 
is partially synchrotron self-absorbed (which leads to a > 0), the 
predi cted co rrelation between the height and the flux was derived 
in lzdziarski] d2012l) . z oc F? Ir , where q = (s + 6)/(2 s + 13) - 1/2. 
Here y is the index of the power-law electron distribution in the 
jet, and F; ntr is the emitted radio flux. As we have shown in Sec¬ 
tion [4] the radio spectrum of Cyg X-3 in the hard state has indeed 
a > 0, and the above formalism approximately applies. Then, at a 
higher flux, and thus a greater height, the photon path through the 
wind becomes shorter, and the orbit-averaged optical depth, f ff , for 
free-free absorption becomes lower. A further effect is an increase 
of the temperature of the wind in response to the increased X-ray 
flux, which also reduces f B . The presence of these effects leads to 
the intrinsic correlation index, /?j ntr , becoming lower than the one 
observed, 


Pmtr = rp. 


d In F[ ntr 
' ~ dlnF R 


1 + 


dfff 

d In F R ’ 


( 2 ) 


where F R is the observed radio flux (ZSPL11). The value of r can 
be estimated based on the dependence of the depth of an orbital 
modul ation of the radio flux on v. In the case of Cyg X-l. lZdziarskil 
j2012l) estimated r = 0.8, leading to p intr ^ 1.4 in that binary. 

In the case of Cyg X-3, the orbital modulation of the hard-state 
radio emission is very weak, with the full depth of the modulation 
of 15 GHz emission of = 5 per cent (Zdziarski et al., in prepara- 
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Figure 6. (a) The colour-colour diagram for selected Galactic BH binaries and Cyg X-3. The soft and hard colours are for the energy flux ratios of the bands 
4-6.4 keV to 3—4 keV and 9.7-16 keV to 6.4-9.7 keV, respectively, (b) The Eddington ratio-hard colour diagram. On both panels, we show Cyg X-l (green 
triangles, 1.86 kpc, 16 Mq), XTE J1650-500 (blue triangles, 2.6 kpc, 7 Mq), LMC X-l (cyan triangles, 50 kpc, 10.9 Mq), XTE J1859+226 (yellow triangles, 
8 kpc, 8 Mq), XTE J1118+480 (magenta triangles, 1.72 kpc, 7.5 Mq), LMC X-l (green squares, 50 kpc, IOMq) GRO J1655—40 (magenta squares, 3.2 kpc, 
6.0 Mq), GX 339—4 (black squares, 8 kpc, 12 Mq), GRS 1915+105 (red open circles, 8.3 kpc, 12.4 Mq) and Cyg X-3 (black stars, 8 kpc, 4 Mq). In parentheses, 
we give the plotted symbols and the adopted distance and mass. 


tion). This is much weaker than that in the hard state of Cyg X-l, 
where it is 30 per cent JZdziarskill2012ll . This weakness is most 
likely due to two effects. First, the separation between the binary 
components in Cyg X-3 is about an order of magnitude lower than 
that of Cyg X-l (given the orbital period of the former is much 
shorter than of the latter). Second, the height along the jet at which 
a given radio frequency is emitted is most likely higher in Cyg X-3 
than in Cyg X-l, due to the accretion rate in the former being much 
higher than in the latter. Therefore, the asymmetry of the paths of 
15-GHz photons through the wind between the superior and infe¬ 
rior conjunction in Cyg X-3 is much lower than in Cyg X-1. This is 
then likely to lead to the value of f ff at 15 GHz being much lower in 
Cyg X-3 than in Cyg X-l, in spite of the wind mass loss rate being 
about an order of magnitude higher. 

Given that, we expect the value of r in Cyg X-3 to be much 
closer to unity than that in Cyg X-1. Thus, the above effect of the in¬ 
trinsic correlation index being lower than the observed one is prob¬ 
ably minor in Cyg X-3. Therefore, it is unlikely that the value of 
p, which we have estimated correcting for the X-ray absorption, is 
significantly higher than our value of =il.5. 

5.2 The nature of the hard state and comparison to other 
X-ray binaries 

The anticorrelation between hard and soft X-rays within the hard 
state of Cyg X-3 is a consequence of the soft spectrum becoming 
softer when brighter. This leads to a pivoting in the =T0-20 keV 
energy range. Then, a flux increase in soft X-rays corresponds to 
a decrease in hard X-rays, and vice versa. This pivoting pattern is 
seen, e.g., in fig. 2 of SZM08, see their spectra 1 and 2. 

On the other hand, Cyg X-l has only moderate spectral vari¬ 
ability in the hard state. Based on the energy depende nce of the 
fractional variability, see, e.g., fig. 12 in lZdziarski et alj ( 120021) . the 


hard-state pivot is at ~ 100 keV. Consequently, there is no observed 
anticorrelation between the soft and hard X-rays in the hard state, 
as found by ZSPL11. 

As calculated in lZdziarski et alj J2002I ). the X-ray variability 
patterns in the hard state are well explained by originating from 
thermal Comptonization. If the dominant driver of the variability 
is in a changing flux of seed photons, the spectrum softens with 
the increase of the total flux, leading to a pivoting behaviour. The 
lower value of the pivot in Cyg X-3 than in Cyg X-l is most likely 
related to the hard state in the former being in fact intermediate. In 
Cyg X-l, the transition from the hard to the intermediate state is 
indeed charact erized by pivoting around ~ 10-20 keV, see fig. 13 in 
IZdziarski et al.i 12002 ). 

Given our results in Section [3] and the discussion in Section 
I5T1 the absorption in Cyg X-3 has to be strong, and the intrinsic 
soft X-ray spectrum in the hard state has to be rather soft, T ~ 2. 
Thus, the hard state of Cyg X-3 corresponds indeed to the inter¬ 
mediate state of BH binaries. As found by ZSPL 11 (see their figs. 
8a, b), the positive correlation between the radio and bolometric 
X-ray fluxes in Cyg X-1 is present in both the hard and intermedi¬ 
ate states, and breaks down only in the soft state. We thus expect 
a similar situation in Cyg X-3, which can explains the observed 
correlation even if its spectral state is intermediate. 

If the apparent hard state in Cyg X-3 is intrinsically intermedi¬ 
ate, this can explain the difference in its spectral shape with respect 
to BH binaries in the hard state without scattering in the stellar 
wind, invoked by lZdziarski et al.l J2010h . The apparent peak at ~10 
keV in the hard-state EF(E) spectrum is then an artefact of strong 
absorption, and the spectrum at higher energies is due to non- 
therm al, rather than thermal, Comptonization jHialmarsdotter et all 
2008). Also, the hard state of Cyg X-3 appears then to be simi¬ 
lar to the hardest spectral state of the BH LMXB GRS 1915+105, 
see, e.g., fig. 4 in lDone, Wardzinski & Gierlinskil^20041 ). where the 
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hardest spectra are similar to the intrinsic hard-state spectra of Cyg 
X-3 shown in Figs.[4jb,d). 

In order to determine the relationship of Cyg X-3 to other 
accreting X-ray binaries, we compare quantitatively their over¬ 
all X-ray and radio properties. We first compare the posi¬ 
tions on the X-ray colour-colour and colour-lu minosity diagrams. 
We u s e the compilations of the data f r om Done & Gierli riskil 
^20031). iDon e, Gierlihs ki & Kubotal d2007t) . iDone et al.l i 20041) and 
Iffladstone et al] J2007I) for Galactic BH and weakly-magnetized 


NS binaries. Those authors used the soft and hard colours defined 
by the energy flux ratio of the bands 4-6.4 keV to 3^4 keV and 
9.7-16 keV to 6.4-9.7 keV, and the bolometric luminosity in the 
Eddington units. We have updated the values of the distances and 
BH ma sses according to t he cur rent best estimates, in particula r 
those of lCasares & Jonkerl d2014h and lCorral-Santana et al.l J20I6I) . 
We calculate the colours and Eddington ratios, Tboi/TEd d, for Cyg 
X-3 assuming <7 = 8 kpc, M = 4Mq dZdziarski et al. ITo 1 3l) and 
using the ASM-BAT absorption-corrected spectra shown in Fig. 
|4jb). Here Le dd is the Eddington luminosity, calculated for pure He 
(X = 0) for Cyg X-3, and for the cosmic composition (X = 0.7) 
for other binaries. To calculate the fluxes within a part of a bin, 
we assume the spectral shape within a bin to be flat in EF(E). The 
spectra at E < 1.3 keV are most likely dominated by disc black- 
body emission, and we thus add a contribution from those energies 
assuming F(E) oc E 113 . We obtain the bolometric fluxes of 2.3, 5.5 
and 8.0 x 10~ 8 erg cnr 2 s _1 in the hard, intermediate and soft state, 
respectively, which correspond to the bolometric (isotropic) lumi¬ 
nosities of 1.8, 4.2, 6.2 x 10 38 (D/8kpc) 2 erg s _1 , and the Eddington 
ratios of 0.18, 0.42 and 0.62(D/8 kpc) 2 (M/4MQ) _1 , respectively. 
We can see in Fig. [6] that indeed the colours and Eddington ratios 
of the spectral states of Cyg X-3 (at 4Mq) are wit hin the range oc¬ 
cupied by GRS 1915+105 (using the data of lDone et alj|2004l . on 
both the colour-colour and colour-Lboi/TEdd diagrams. 

On the other hand, iHialmarsdotter et ail d2009l) used similar 
diagrams and argued that the most likely BH mass is 30Mq. In 
fact, that mass may also be compatible with our values of colour- 
Lb 0 i/LEdd; an d it would correspond to moving the points on Fig. 
[6jb) down to a region populated by other BH binaries, but not GRS 
1915+105. However, the Cyg X-3 points on the colour-colour dia¬ 
gram would still be within those of GRS 1915+105, which repre¬ 
sents an inconsistency. Als o, suc h a high B H mass is not found 
in other BH binaries (see [Lavcock. Maccarone & Christodouloul 
12015l for the lack of dynamical evidence for a high mass BH in 
IC10 X-l), as well as it is in contradic tion with the later results o f 
IZdziarski et alj d2013l) . Compared to IHialmarsdotter et al.l J 20091) . 
the absorber in Cyg X-3 adopted by us is stronger, which would 
move the values of the colours in that paper to lower values, making 
them more compatible with those of GRS 1915+105. We also note 
that the distance to GRS 1915+10 5 ha s been revised down from 
12.5 kpc to 8.6^ 2 ® kpc dReid et al .1120141) . which has decreased its 
Eddington ratios by a factor of 2.1 (at the best-fit distance). This 
results in our Cyg X-3 points for M = 4Mq to be entirely within 
the region occupied by GRS 1915+105. 

We have also compared the colours and the Eddington ratios 
of Cyg X-3 to those of weakly magnetized NS binaries, assum¬ 
ing now M = 1.5Mq. We have found the results of this compar¬ 
ison compatible with Cyg X-3 containing an NS, with its colours 
and Eddington ratios similar to those of Z sources. However, the 
strength of the radio emission relative to the X-ray in Cyg X-3 ex¬ 
ceeds that of all of BH bin aries, see below, and it i s much much 
higher than for NS ones (Migliari & Fended 1 20061 : ICorbel et al.l 
l2013al) . 



Figure 7. The correlation between the unabsorbed bolometric luminosity 
in Eddington units with the radio 15 GHz luminosity for Cyg X-3 (8 kpc, 
4Mq; black error bar's) compared to that of Cyg X-l (1.86 kpc, 16Mq; 
the red triangles, green squares and blue circles for the hard, intermediate 
and soft state, respectively, from ZSPL11), GX 33SM- (8 kpc, 12Mq; green 
squares fitted by the green dashed line, from lZdziarski et al .2(11 but mul¬ 
tiplied by 1.25 to account for those measurements being at 8.6 GHz) and 
GRS 1915+105 (8.6 kpc, 12.4Mq; cyan squares with error bars), with the 
magenta squares corresponding to the plateau \ state using the criteria de¬ 
scribed in Scction [5.2l 


We then compare the radio vs. the bolometric luminosity re¬ 
lationship of Cyg X-3 with those of selected BH binaries, see Fig. 
[7] where the y axis shows vL v /Leh at 15 GHz. The black error 
bars shows the points for Cyg X-3 at M = 4Mq and <7 = 8 kpc, 
where the shown values of L bol have been corrected with respect to 
L(l.3-100keV) in the same way as described above for Fig.[6jb). 

We compare first Cyg X-3 with GX 339- 4, for w hich the 
radio/X-ray correlation is the best studied one ICorbel et alj2013ah . 
Unfortunately, the published correlations are given for narrow 
bands of either 3-9 keV or 1-10 keV only. As we have argued 
above, a narrow-band X-ray flux is likely to be significantly af¬ 
fected by spectral variability (as it is very significantly the case 
in Cyg X-3, and significantly in Cyg X-l, see ZSPL1T). The only 
radio/Lboi correlation available for an LMXB appears to be that 
of IZdziarski et alj 12004!) for GX 339^4, which we show in green 
points. The magenta line shows the best fit power-law dependence, 
with p = 0.79 ± 0.07. The radio fluxes are at 8.6 GHz l lCorbel et al.l 
1 200 3), which we multiply by a factor of 1.25 in order to estimate 
the 15 GHz flux using the average a = 0.4 observed in the hard 
state of GX 339^4 dCorhel et alj|20 1 3bl) . On the other hand, we 
have checked that the 8.3-15 GHz index of Cyg X-3 is = 0 in its 
hard state , cf. Figs. [5ta-b), and it is — 0 for Cyg X-l in the hard 
state iFender et alj|2000l) . Thus, the effect of the radio band on the 
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presented correlation appears minor. We see Cyg X-3 in its hard 
state is significantly more radio-loud than GX 339^4, by a factor of 
—4 on average (though there is a systematic uncertainty concerning 
the correction for the intrinsic X-ray absorption in Cyg X-3). 

As we have shown above, Cyg X-3 appears relatively similar 
spectrally in X-rays to GRS 1915+105. In order to them in the ra¬ 
dio band, we use the intense Ryle telescope monitoring of GRS 
1915+105 conducted during MJD 49856-53898, and the ASM 
data, during MJD 50088-55859. Given the fast variability of that 
object, we have tested the effect of various averaging methods. We 
have finally chosen to average the radio over observation periods 
containing gaps shorter than 0.2 d, which resulted in the average 
length of an averaged segment of 0.09 d (usually one per day) 
with a — 0.07 d. We then correlated those data with the ASM 
dwell data which were within ±0.2 d from the average time of a 
given radio segment. We have corrected the averaged ASM data for 
line-of-sight absorption assum ing Nh = 3.8 x 10 22 cm 2 febisawal 
Il998l : lchapuis & Corbelll2004l) and the intrinsic spectra with T = 2 
within each ASM bin. We have further added the contribution from 
E < 1.3 keV in the same way as for Cyg X-3, i.e., assuming 
F(E) oc E 1 ' 3 , see above. We have also considered the available 
BATSE data for this objecJJ We have found that addition of the 
20-160 keV fluxes to those of the ASM was possible only for about 
40 per cent of the ASM points, and an upper part of the correlation 
was not covered by the ASM+BATSE points. The average 20-160 
keV to 1.3-12.2 keV flux ratio equals 0.11. Given that the EF(E) 
spectra of GRS 1915 +105 in al most all states are decreasing above 
several keV l lDone et alj|2004l ), we estimate the flux contribution 
from 12-20 keV to be similar to that from > 20 keV. We have thus 
applied an approximate bolometric correction of a factor of 1.2 to 
the fluxes from the ASM. 

We show the results by the cyan error bars in Fig. [7] We see 
GRS 1915+105 doe s not show a distin ct radio/X-ray correlation. 
On the other hand, iRushton et alj 1201(1) claimed a strong correla¬ 
tion within the plateau (steady and hard) state of GRS 191 5+105. 
When using the same time intervals as those given in lRushton et al.1 
feOld) as used for their correlation, we have found a rather scat¬ 
tered diagram, with a steep positive slope at high radio fluxes and a 
negative one at low radio fluxes. When we used the same criteria as 
given in that paper, i.e., the ratio of the count rates in the ASM third 
to second channels of > 1.5 and the total ASM rate of 30-50 s -1 , 
we found a similar discrepancy. We have recovered a correlation 
only after imposing an additional criterion of the radio flux of > 40 
mjy (which, however, implies that the X-ray state by itself does not 
set the radio behaviour, in contrast to Cyg X-3). Those points are 
shown in Fig.[7]by green crosses. We do see some correlation, with 
the slope of 2.1 ± 0.2, somewhat steeper than 1.7 ± 0.3 given in 
IRushton et al.l d20ld) . 

We see that GRS 1915+105 occupies a similar range of 
F-boi/TEdd as Cyg X-3 (as also seen in Fig. [6jn), but it is both less 
radio loud and with a much less clear radio/X-ray correlation. The 
likely reason for that is the difference in the jet environments. Cyg 
X-3 is a binary with a high-mass donor, and interaction between its 
wind and the jet le ads to presence of shocks^ in pa rticular of rec- 
ollim ation shocks fYoon & Heinz 201 5|; lYoon, Zdziarski & Heinzl 
l2016h . Those shocks are likely to efficiently accelerate electrons, 
leading to enhanced radio emission. However, no such interaction 
takes place in GRS 1915+105, an LXMB. 

On the other hand, Cyg X-l, which is also a high-mass X-ray 


3 http: //gammaray.nsstc.nasa. gov/batse/occultation/ 


binary, has the hard-state radio fluxes, shown by the red triangles in 
Fig. [7] an order of magnitude below those of the GX 339^4 correla¬ 
tion. An explanation for that may be provided by the wind ram pres¬ 
sure around the jet in Cyg X-l being much lower than that in Cyg 
X-3, owing to both the lower mass-loss rate and the larger separa¬ 
tion. The jet lateral pressure in Cyg X-l may then be always larger 
than th at of the wind, provided the jet power is above some critical 
value jYoon et alj 20161) . Then no recollimation shock would form. 

Furthermore, as discussed in Section [5Tl the radio flux in Cyg 
X-l is significantly free-free absorbed in the wind from its donor, 
see ZSPL11 and IZdziarskil J2012I ). In the model of the latter, the 
slope of the intrinsic hard-state correlation in Cyg X-l is close to 
that of the lower branch of the radio/X-ray correlation, first found 
for HI743-322 fcoriat et al.l201lh . However, the correction for the 
wind absorption in Cyg X-l is relatively uncertain. More detailed 
knowledge of the wind structure in Cyg X-l is needed to determine 
the degree of free-free absorption. 

5.3 Consequences for radio/X-ray correlations in other 
objects 

Our results point out to the importance of taking into account an 
estimate of the bolometric flux when correlating X-rays with radio, 
reinforcing the similar conclusion of ZSPL11. Due to the pivoting 
spectral variability in X-rays the hard X-rays are often a very bad 
predictor of the bolometric flux, and can be either uncorrelated or 
negatively correlated with it. This effect appears to be present also 
in Seyfert galaxies, in which a common variability pattern in soft 
X-rays is the steepening of their soft X-ray spectrum with the in¬ 
creasing flux. This then leads to a pivot in hard X-rays. 

This effect can ex plain th e discrepancy between the result of 
ISalvato, Greiner & Kuhlbrodl J2004l) , who found the radio and soft 
X-ray emission of thei r sampl e of Seyferts to be strongly correlated, 
and the result of Burlon et alj l l2013l) . who found a lack of such cor¬ 
relation for hard X-rays, using SwiftfBAT data. The latter authors 
doubted then the reality of the soft X-ray correlation. However, this 
discrepancy can be resolved by a pivoting spectral variability of X- 
rays. In Cyg X-3, it leads to the 10-20 keV flux to be completely 
uncorrelated with radio, see Fig.[2jg), in spite of the correlation in 
soft X-rays, Fig.[2je), being very strong. 


6 CONCLUSIONS 

We have studied implications of the presence of the positive and 
negative correlation with the radio flux for the soft and hard X- 
ray fluxes, respectively, in the hard spectral state of Cyg X-3. 
We have found that the radio flux is then completely uncorrelated 
with the observed, absorbed, broad-band, 1.3-100 keV X-ray flux, 
which appears to be a good approximation to the absorbed bolo¬ 
metric flux. An intrinsic lack of correlation is unlikely on theoreti¬ 
cal grounds, as well as it is unlikely given the presence of the strong 
narrow X-ray band correlations. 

However, we can recover a clear positive correlation if the 
absorption by the stellar wind of the donor in Cyg X-3 is strong 
enough. This implies that the intrinsic soft X-ray spectrum in the 
hard state is relatively soft, similar to those of intermediate states 
of BH binaries. This nature of the apparent hard state of Cyg X-3 
is also consistent with the luminosity of the hard state being much 
higher in this object than in the hard state of the high-mass BH 
binary, Cyg X-l. 

We have found the radio spectra in the hard state to be hard on 
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average, with a 0. Such spectra are characteristic to BH X-ray 
binaries in the hard state. 

We have compared Cyg X-3 with other X-ray binaries. On the 
X-ray colour-colour and colour-Eddington ratio diagrams, Cyg X- 
3, after correcting its spectra for absorption and for the BH mass 
of ~ 4Mq, appears similar to GRS 1915+105, with similar colours 
and Lboi/^Edd ~ 0.1-1. On the other hand, Cyg X-3, at 1.5Mq is 
also similar on such diagrams to NS Z-type binaries. However, the 
radio emission of those sources is much weaker than that of Cyg X- 
3, which is both the most radio-luminous and, in its hard state, the 
most radio-loud (i.e., with the highest ratio of the radio to bolomet- 
ric luminosities) source among all known X-ray binaries. Still, its 
radio loudness is much closer to those of BH binaries than to those 
of NS ones. This represents a strong argument for the presence of a 
BH in Cyg X-3. We interpret the radio loudness and radio fluxes of 
Cyg X-3 being much higher than those of GRS 1915+105 as due to 
an enhancement of the jet emission in the former due to interaction 
of the jet with the very strong stellar wind from the donor. 

We have found (Appendix 0 that the flux distributions of the 
radio emission and soft X-rays can be described by two log-normal 
functions, one corresponding to the hard state and one to soft state. 
This indicates a multiplicative and correlated character of variabil¬ 
ity of Cyg X-3 in both the hard and soft state and in both the accre¬ 
tion flow and the jet. 

Finally, we have presented (Appendix 0 the long-time light 
curves of Cyg X-3 using the available data, and compared them 
to the occurrences of high-energy y-ray detection. We have found 
that a low hard X-ray flux, F/(F) S; 0.1, is both a necessary and 
sufficient criterion for detectable y-rays. We have found that this 
criterion has not been satisfied since MJD 55650 (2011 March 30) 
until now, MJD 57333. 
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APPENDIX A: THE FLUX DISTRIBUTIONS AND LIGHT 
CURVES OF CYG X-3 

We show here distributions of the fluxes and the light curves for the 
X-ray monitors used here and for the Ryle/AMI. Fig JAll shows the 
histograms of the 15 GHz, ASM, MAXI and BAT fluxes and count 
rates (for simplicity, we use the s ymbol F for both). The histo gram 
for the BATSE flux is shown in lHialmarsdotter et al.l 120081) . The 
rates for 15 GHz, and ASM and MAXI are fair estimators of the 


jet power in the former case and the accretion power in the latter, 
see Section[3] On the other hand, the hard X-rays are anticorrelated 
with the accretion luminosity in the hard/intermediate state due to 
the spectral pivoting. Thus, their distribution does not reflect the 
accretion power. 

Consequently, we find the 15 GHz, ASM and MAXI his¬ 
tograms can be fitted by a sum of two log-normal distributions, 
one present in the hard/intermediate state and one in the soft state, 
while it is not the case for the BAT histogram. The values of 
the mean and standard deviation for the hard/intermediate-state 
and soft-state log-normal distributions for the 15 GHz, ASM and 
MAXI data are (-2.00,0.22; -2.09,0.90), (0.84,0.13; 1.28,0.11), 
(-0.40,0.06; -0.01,0.04), respectively (in log 10 units). We see that 
the width of the flux distribution in each of the state in the ASM 
data is twice as large as that of the MAXI data. This is due to the 
MAXI 10-20 keV channel containing the pivot of the spectral vari¬ 
ability, see Fig.[2jg), which reduces the range of the flux variability. 
On the other hand, the ASM 5-12.2 keV channel is still below the 
pivot, see Fig. |2jc). All of the ASM and MAXI distributions are 
narrower than the corresponding 15-GHz ones, especially in the 
soft state. This reflects the dramatic variability of the radio emis¬ 
sion observed in that state. 

The presence of a log-normal distribution shows that the un¬ 
derlying physical process has to be multiplicative, with variations 
coupled together on all time scales l lUttlev. McHardy & Vaughan! 
2005 1. Interestingly, we find the log-normal variability to be present 
in the both main sta tes of Cyg X-3, hard and soft. On the other hand , 
lUttlev et .ill ( 2005 1 and IPoutanen, Zdziarski & Ibragimov! ( 2008 ) 
have found it in the hard state of Cyg X-l on ~0.1-10 s (for soft 
X-rays) and ;S 0.1 d (for soft X-rays and radio) time scales, respec¬ 
tively. However, it appears that a log-normal flux distribution has 
not been found before in the soft state of an X-ray binary. 

Fig. lA2l shows the light curves. Since the variability in all en¬ 
ergy bands spans a large ratio and it is log-normal for some de¬ 
tectors, we normalize each light curve to its geometric unweighted 
average over the respective observation time span, i.e., we calcu¬ 
late (In F). For that, we have to select only points with the measured 
fluxes of Fj > 0. This will lead to a certain overestimate of the over¬ 
all average, but the effect is minor for the light curves presented be¬ 
low. For notational simplicity, we denote (F) = exp (In F ). Also, in 
order to reduce a contribution from the measurements of poor qual¬ 
ity, we include only data with the fractional error of AF,/F, < 0.5, 
where AF, is the flux error. We show the 1-d average rates except 
for the BATSE, where we show the 3-d average. In the case of BAT, 
we use the public 15-50 keV dattQ (Krimm et ail2013l) . That light 
curve is very similar to our 14-50 keV one. For the MAXI, we show 
the 2-10 keV light curve instead of the full 2-20 keV one in order 
to have as much correspondence as possible to the ASM 1.3-12.2 
keV light curve. 

The average rates and standard deviations are given in Table 
Al. The standard deviations correspond to the excess variance, i.e., 
the variance due to the measurement errors is subtracted, 

0-2 = x 2 K Fi ~ <ln F>2 - ( AF 'l F >r\ ■ (ad 

i 

The lower value of <x for the BATSE data than that for the BAT 
is probably due to the former data being averaged over 3-d rather 
than 1-d intervals. The large values of a for all the radio data is the 
consequence of the dramatic radio variability in this object, show- 

4 http://swift.gsfc.nasa.gov/results/transients/CygX-3/ 
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Figure Al. The histograms of the 15 GHz fluxes and the ASM, MAXI and BAT count rates, from top left to bottom right. The first three histograms are well 
represented by a sum of two log-normal distributions, shown by the black solid and red dashed curves, which correspond to the hard and soft state, respectively. 


ing both very bright flares and very weak states, see Fig. IA2I For 
the 15 GHz data, we have used the estimate of A Fj/Fj = 0.1 (Sec- 
tion[2]l. For clarity, we do not show the GBI light curves, which are 
relatively similar to the Ryle/AMI one. 

We also mark the intervals corresponding to the de¬ 
tected high-energy y-ray emission in Fig. IA2I These are 


MJD 54566-54647, 54821-54850, 55324-55326 i Bulgarelli et al 
20121). 54750-54820, 54990-55045 dFermi LAT Collaboration 


2009 1. 55343-55345 ^Williams etaTll2011 


IBulgarelli et alj|2012h . 


ulg; 

55586 -55610, 55642-55644 dCorbelet al. 1)20 121) and MJD 57398- 
57412 jLoh et af1l20 1 6l : iTavani et alJl2011 ). All those periods cor¬ 
respond to the soft state as defined in our paper, and, with exception 
for the last occurence, to periods with both the lowest 15-50 keV 
flux detected, with F/(F ) S; 0.1, and occurrences of radio flares, 
wi th the 1-d averaged 15 GHz flux of F/(F ) > 5 (as noted, e.g., 
bv ICorbel et al.ll2012l : IBulgarelli et al.ll2012h . However, the former 
of these conditions, i.e., a low hard X-ray flux, appears to be both 
necessary and sufficient. We note that this condition has not been 
satisfied during MJD 55650-57397, see Fig. IA2I Indeed, no high- 
energy y-ray detection from Cyg X-3 has been reported during that 
period. 


Table Al. The logarithmically-averaged detector rates and standard devia¬ 
tions for Cyg X-3. 


Detector 

Energy range 

Unit 

(F) 

(T 

GBI 

2.25 GHz 

mJy 

99 

0.91 

GBI 

8.3 GHz 

mJy 

125 

0.84 

Ryle, AMI 

15 GHz 

mJy 

103 

0.94 

RXTE/ ASM 

1.3-12.2 keV 

count s -1 

10.5 

0.58 

/SS/MAXI 

2-10 keV 

photon cm -2 s -1 

0.350 

0.50 

Swift /BAT 

15-50 keV 

count cm -2 s -1 

0.0336 

0.54 

CGRO/BATSE 

20-100 keV 

keV cm -2 s -1 

1.38 

0.30 
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Figure A2. The long-term light curves of Cyg X-3 normalized to the respective average over the observation length. The cyan, blue, red, magenta and black 
symbols show the rates normalized to the respective average for the ASM (1.3-12.2 keV), MAXI (2-10 keV), BAT (15-50 keV), BATSE (20-100 keV) and 
Ryle/AMI (15 GHz). For clarity of display, the error bars are not plotted. The dashed green line corresponds to the average. The heavy magenta horizontal 
lines correspond to the detected 8 occurrences of high-energy y-ray emission. 
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Figure A2. Continued. 
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Figure A2. Continued. 
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